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D ■ Abstract. 



The k = 3 Read-Rezayi (RR) parafermion quantum Hall state hosts non- 
Abelian excitations which provide a platform for the universal topological quantum 
computation. Although the RR state may be realized at the filling factor v = 12/5 in 
GaAs-based two-dimensional electron systems, the corresponding quantum Hall state is 

■ weak and at present nearly impossible to study experimentally. Here we argue that the 

RR state can alternatively be realized in a class of chiral materials with massless and 
massive Dirac-like band structure. This family of materials encompasses monolayer 
and bilayer graphene, as well as topological insulators. We show that, compared to 
GaAs, these systems provide several important advantages in realizing and studying the 

■ RR state. Most importantly, the effective interactions can be tuned in situ by varying 

the external magnetic field, and by designing the dielectric environment of the sample. 

| This tunability enables the realization of RR state with controllable energy gaps in 

• different Landau levels. It also allows one to probe the quantum phase transitions to 

other compressible and incompressible phases. 
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1. Introduction 

Strongly correlated phases of electrons confined to move in the plane, subject to a 
perpendicular magnetic field, have attracted significant attention since the discovery 
of fractionally quantized Hall conductivity pp. The profound role of topology in this 
extreme quantum limit leads to the presence of quasiparticles that carry a fraction of 
electron charge [2], and fractional (Abelian, or possibly non-Abelian) statistics [21 [3]. 
The prospect of excitations possessing non-Abelian statistics has motivated different 
schemes for topological quantum computation [4\ based on these systems. 

These remarkable phenomena occur in the fractional quantum Hall (FQH) regime, 
when the number of electrons, N e , is comparable to the number of magnetic flux 
quanta N$ through the two-dimension electron system (2DES). Correlated FQH liquid 
states appear at certain partial filling v = N e /Nq> of the active Landau level (LL). 
In traditional semiconductor heterostructures, the physics of a partially-filled n — 1 
LL differs significantly from that of n = LL, due to the node in the single-particle 
wavefunction [5]. As a consequence, the hierarchy /composite fermion states [61 [7J, 
ubiquitous in the lowest Landau level (LLL), are significantly weakened in n = 1 LL, and 
some of the more exotic states, such as the Read-Rezayi (RR) parafermion states [8], 
are likely to be favored. A number of studies have focused on the simplest, k = 2 
non-Abelian member of the RR sequence - the Moore- Read (MR) "Pfaffian" state [3]. 
The Pfaffian, or perhaps its particle-hole conjugate version - the anti-Pfaffian [9], with 
the same topological properties, are widely believed to describe the FQH plateau at 
v — 5/2 [10]. Quasiparticles of the MR and higher RR states obey the non-Abelian 
statistics [3] which is of interest for topological quantum computation jl]. However, for 
the purpose of universal topological quantum computation, MR state is not sufficient, 
and one must go to a higher, k = 3 member of the RR sequence (see Ref. jl] and for 
a general case Ref. [H]). In GaAs systems, experiments [121 [13] have detected only a 
weak plateau at v = 12/5, which theoretical work [8] has tentatively identified with 
the particle-hole conjugate of the k = 3 RR state. This state, however, is very fragile 
and has been seen only in a small number of samples, triggering alternative theoretical 
proposals for its origin [14]. Numerical calculations suggest that the appearance of the 
RR k = 3 state is linked to the finite- width of the 2DES in GaAs [151 HE], but little is 
known about the stability of k = 3 RR state as the interaction is tuned away from the 
bare Coulomb point. For example, in case of k = 2 (Moore-Read) state, it is possible 
to construct a 2-body interaction [17], resulting from particle- hole symmetrization of 
the hard-core 3-body repulsion, which yields a ground-state with high overlap with the 
exact ground state of Coulomb interaction. This model interaction consists only of V\ 
and V3 Haldane pseudopotentials (to be defined below). In case of k = 3 RR state, 
such an approximate 2-body hard-core interaction has not been constructed, but one 
could expect that it requires higher order Haldane pseudopotentials (and perhaps some 
3-body terms as well). 

One of the main disdvantanges of the GaAs-based devices is that their 2DES is 
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buried inside a larger, three-dimensional structure. This unfortunate fact fixes the 
effective interactions at values that are often not optimal for some of the most interesting 
FQH states, including the RR series. For example, the MR state is found to lie very- 
close to the boundary with a compressible phase [TU 19J. Another problem stems from 
the strong dielectric screening and finite well-width [20] in GaAs, which weaken the 
electron-electron interactions and make FQH states fragile. This has been a major 
obstacle in the studies of the non-Abelian states, which could only be observed in ultra- 
high- mobility samples [10J. Thus, it is desirable to find an alternative high-mobility 
2DES with strong effective Coulomb interactions that are adjustable in a broad range. 

Recently, a new class of such high-mobility 2DES which host chiral excitations 
with non-trivial Berry phases, has been discovered. These "chiral" materials include 
graphene and bilayer graphene [21], and, more recently, topological insulators [22], as 
well as certain quantum wells [23]. The chiral nature of the quasiparticles gives rise 
to new electronic properties, including an unusual LL sequence, anomalous Hall effect, 
and suppression of weak localization [21] , When these chiral materials are subject to a 
perpendicular magnetic field, the kinetic energy quenches into discrete Landau levels, 
similar to the usual semiconductors with non-chiral carriers. However, due to the chiral 
band structure |24j and the fact that the surface of these materials is exposed [25], they 
offer new possibilities to tune the effective interactions and explore strongly correlated 
phases. 

In this paper we analyze two practical ways of tuning the effective interactions in 
chiral 2DES, and the effect this has on the non-Abelian FQH states. As a case study, 
we choose to focus on the k = 3 RR state. In a stark contrast to GaAs, we find that the 
chiral 2DESs allow for a more robust RR state along with the possibility of its realization 
in several LLs. Additional insights can be obtained by driving transitions between the 
RR state and the Abelian hierarchy state, or the compressible stripe phase [26]. Such 
transitions can be implemented in chiral and massive 2DES by varying the external 
field. Overall, the tunability of the chiral 2DES allows one to explore a larger region of 
the effective interactions than has been achieved in GaAs. 

2. Model 

We consider a family of 2D materials that are characterized by non-trivial Berry phases. 
One such material is monolayer graphene, a high-mobility atomically thick 2DES |21j, 
where recently several FQH states of the type v — m/3 have been discovered [27J. A 
closely related material, bilayer graphene [21] has similarly high mobility, and exhibits 
interaction-induced quantum Hall states at integer filling factors at low magnetic 
fields [31]. Graphene and its bilayer are characterized by Berry phase 7r |28j (graphene- 
like) and 2-7T [29] (bilayer-graphene-like with an energy gap), respectively. 

Much of previous theoretical work on chiral materials has been restricted to 
graphene, exploring the consequences of the four-fold LL degeneracy (valley and spin) 
that leads to new SU(2) and SU(4)-symmetric fractional states [30J. Instead, we focus 
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Figure 1. (Color online) Ratio V1/V3 of the Haldane pseudopotentials as a function 
of system's parameters: the width w in case of the infinite quantum well model (a), 
the band structure parameter 9 (b), and screening parameter a and distance of the 
dielectric plate d (c). Chiral 2DES offer the possibility to tune V1/V3 in a much broader 
range than it is possible within a quantum well model. 



on the high- field limit, neglecting the multicomponent degrees of freedom, and examine 
the effects originating from the interplay of the Coulomb interaction and band structure. 
We consider a family of band structures introduced in Ref. [24J, which describe a 
number of high-mobility materials, including graphene with massive carriers (mass is 
generated either spontaneously, or as a result of sublattice symmetry breaking [2~T]). 
topological insulators |22j, bilayer-graphene [21] . trilayer graphene [31], and similar 
materials. Pristine graphene, which hosts massless Dirac-like fermions, is contained 
in this model particular case. 

More explicitly, we consider a family of 2 x 2 Hamiltonians with Berry phase n and 
27t [24] . The case of 7r-carriers is realized in graphene, topological insulators, special 
quantum wells |23j ; the case of 27r-carriers occurs in bilayer graphene, where the energy 
gap can be controlled in a wide range by a perpendicular electric field [21] . These 2x2 
effective Hamiltonians can be derived pertubatively starting from a tight-binding model 
for each of the materials with a given chirality A, however the effective model will be 
applicable in a narrower energy range for larger A [32]. As it is well known for the case 
of graphene, the one-body eigenstates of such Hamiltonians possess a spinor structure. 
Consequently, the harmonic oscillator raising/lowering operators will mix the different 
components to produce an effective form-factor that has a richer structure than in case 
of a simple parabolic band such as in GaAs, making it possible to tune the effective 
interaction. 

For 7r carriers, the single-particle wavefunctions are given by [23] ip n = 
(cos8 n <f)\ n \_i,sm8 n <j)\ n \), where <fr n is the wave function of the nth non-relativistic LL 
(n = ±1, ±2, . . .), and parameter 9 depends on the ratio A/(^o/^s) ) where A is the 
mass gap, v the Fermi velocity, and £b = ^/hc/eB the magnetic length. We use the 
notation (crir, n) to denote the nth LL for air carriers (a = 1, 2). As a consequence of the 
spinor wavefunction, the effective form factor [5] F£(q) that describes the interaction 
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projected to a (7r, n) LL is given by 

FM = cos 2 6L lnhl (t3>>) + sin 2 9L H , (1) 

where is the fcth Laguerre polynomial, and for simplicity we omitted the index of 
9. The form-factor is a mixture of the (|n| — l)th and |n|th LL form-factors in a non- 
relativistic 2DES with parabolic dispersion. At 9 — 7r/4, the above equation reduces to 
the form-factor of graphene [33], however, quite generally by varying A/(Hvq/£b), one 
can realize any value of 9 G (0; vr/2). 

Similarly, for carriers with Berry phase 2ir, the single-particle wavefunctions are 
ip n = (cos0 n 0| n |_i, sin^^ini+i), and the form-factor is a mixture of standard (\n\ — l)th 
and (| rz. j + l)th form-factors, 

F?(q) = cos 2 9L\ n \_i (f^j + sin 2 9L H+l . (2) 

The tunable form of the effective interactions Eqs.( 1Tf2"j) provides a way to engineer 
transitions between strongly correlated phases in situ by changing the field. As pointed 
out by Haldane [5], an interacting many-body system of electrons, confined to a partially 
filled LL, is defined by a finite set of numbers - the Haldane pseudopotentials. For 
rotationally invariant systems, these numbers V m represent the amplitudes of a state 
of two electrons with the relative angular momentum m. The physics of the FQHE is 
determined by small-m pseudopotentials, i.e. V±, V 3 ,V 5 , . . .. The role of higher V^'s is 
to enforce charge neutrality, but they do not affect the incompressibility. Therefore, 
the ratio of the two strongest pseudopotentials, V1/V3, can be a rough indicator of the 
expected many-body phases. In Fig. [H(b) we show this ratio for several types of charge 
carriers and LLs. The variation of V1/V3 as a function of 9 is to be contrasted with 
FigH^a) where the same quantity is plotted for an infinite GaAs quantum well of width 
w. In GaAs the width of the well can be tuned by electrostatic gates; however, it 
invariably reduces V1/V3, whereas the chiral 2DES allow for both increase and decrease 
of V1/V3, and in a wider range. 

Up to this point, we have discussed the change in the effective interaction 
V(q)\F(q)\ 2 that resulted from a modification of the single-particle wavefunctions, and 
therefore F(q). The second method, proposed in Ref. [25], allows one to directly change 
V(q). Chiral materials, such as graphene, are often exposed to the environment, which 
allows for dielectric material to be deposited on top of them. We consider a setup where 
graphene sample is situated in a dielectric medium with permittivity ei, and a semi- 
infinite dielectric plate with permittivity e 2 7^ ei is placed at a distance d/2 away from 
the graphene sheet. The effective interactions between electrons in graphene change due 
to the surface charges induced at the boundary between dielectrics [25J: 

V(r) = h a , where a = — -. (3) 

The ratio djls controls the effective interactions within a partially filled LL. However, 
the overall energy scale is also modified and this has an impact on the magnitude of 
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the excitation gap. The gap should be multiplied by a factor eo a As/ e i if comparison is 
to be made with GaAs 2DES. Again, an important advantage of this setup is that the 
interactions can be tuned in situ by varying the magnetic field B, which modifies the 
ratio djls- The consequences for the many-body system are illustrated in Fig. QJc). In 
the dielectric case, V1/V3 is found to have a substantially larger variation as a function 
of a than the quantum well model or the tuning of the band structure through 9. In 
particular, for intermediate d ~ £3 and negative a, V1/V3 can be readily increased by 
50% compared to its bare value (a = 0). 

Therefore, based on the behavior of V1/V3, we expect that chiral 2DES offer many 
possibilities of tuning the interactions. By increasing this ratio, we expect to stabilize 
the incompressible phases; on the other hand, by reducing it, one can drive a transition 
to compressible phases, such as stripes and bubbles. In the following, we will explicitly 
verify these statements in exact diagonalization calculations, and determine the stability 
of the k = 3 RR state. In the operational sense, one can define several criteria for the 
"stability" of a state; the more of those criteria are met, the more "stable" the state 
is. For example, one can compare the overlap (scalar product) between an exact many- 
body ground state and the RR wavefunction; if the overlap, defined below, is consistenly 
close to unity for a number of system sizes, we consider it as one of the indicators that 
the state is in the RR phase. However, there are known examples of different phases of 
matter which possess ground states of high overlap with each other in finite systems [35J ; 
conversely, a relatively poor overlap cannot be taken as a definitive proof that a phase 
is not realized. Additional indicators to identify the phase of matter are, for example, 
the ground state degeneracy on the torus, and the shift on the sphere (selected to 
minimize the total ground-state energy). Finally, the robustness of a phase is also 
determined by the magnitude of the gap for creating neutral and charged excitations. 
All of these indicators should be assessed together in determining the nature of a phase 
being studied. 

3. Method and results 

We consider an interacting, iV-electron problem with a Landau-level spectrum given 
above, on a compact surface threaded by N$ magnetic flux quanta using exact 
diagonalization. In the study of FQHE, two kind of surfaces are available that preserve 
the translational invariance of an infinite 2DEG: sphere [6] and torus [36]. The two 
choices of boundary conditions illustrate the specific features of a FQH state under 
investigation: on a sphere, the FQH state couples to the curvature of the manifold, which 
is characterized by the topological number called shift [37]. The shift produces a small 
offset between N$> and the magnetic monopole, whose strength is denoted by 2S, placed 
in the center of a sphere. As a consequence, different FQH states "live" in different 
Hilbert spaces and in principle can be directly compared only after extrapolation to the 
thermodynamic limit. On the other hand, the flat surface of a torus leads to a unique 
definition of iV$ for given N and u, so that different candidate states describing the 
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Figure 2. (Color online) Overlap between the k — 3 RR state and the exact ground 
state in (tt, 1) and (tt, 2) LLs as a function of 9, for systems of different sizes N. High 
overlap is found close to the bare Coulomb interaction in n = 1 non-relativistic LL. 
At cos 2 9 0.15, the ground state undergoes a change in shift, and a hierarchy state 
sets in (right). For cos 2 9 < 0.8 in (tt, 2) LL, the system undergoes an incompressible- 
compressible transition to a stripe phase. 



same filling u are all realized in the same Hilbert space. The caveats of this geometry 
are the additional geometric parameters, the angle and aspect ratio of the torus; the 
Hamiltonian of a finite-size system depends on these parameters and their specific values 
favor one FQH phase over the others. Thus, an analysis is slightly more involved but 
the advantage is that a ground-state degeneracy can be used to identify the phase. For 
example, in order to detect a stripe phase, we use a rectangular geometry with a specific 
aspect ratio that best accommodates the stripe in one direction. In contrast, on the 
sphere FQH states always appear as non-degenerate, zero angular momentum ground 
states. 

We will mostly present data for the spherical geometry because this allows for a 
number of simple diagnostic quantities to be evaluated in a straightforward manner. 
The spherical representation of broken-symmetry states, however, necessarily contains 
defects, and periodic boundary conditions have to be used in this case. A number of 
useful insights can be inferred from the study of the energy spectrum; in addition, we 
will use overlap calculations to compare an exact, many-body ground state ^ exact, with a 
numerical representation of a trial wavefunction, Atrial- The overlap is defined as a scalar 
product between two normalized vectors, O = | (^triail^cxact) I- If O is consistently close 
to unity for a number of system sizes considered, we consider the trial wavefunction 
to be a faithful description of a FQH phase. From the knowledge of a ground-state 
wavefunction \I/ , we also evaluate the (projected) static structure factor [39J. Sharp 
peaks in the structure factor indicate the onset of compressible phases [40J. 

In Fig. [2] we show the overlap between the k = 3 RR wavefunction and the exact 
ground state of a system in (tt, 1) and (tt, 2) LLs, as a function of 9. We set the shift 
is —3, corresponding to that of the RR state, and do not consider other candidate 
states with different values of the shift [38]. RR model wavefunction can be obtained by 
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Figure 3. (Color online) Overlap between the k — 3 RR state and the exact ground 
state in (2tt, 1) LL as a function of 9, for various system sizes N. The region of high 
overlap is found in the window between the stripe phase and the hierarchy phase. 



diagonalizing a four-body short-range interaction [8] , which quickly becomes numerically 
prohibitive as the number of particles grows, or recursively via Jack polynomials jH]. 
The filling factor is fixed at v = 3/5, but in the absence of LL mixing, our results 
directly apply to v = 2/5 via particle-hole symmetry The highest overlap is found for 
the effective interaction that is in the vicinity of the pure n = 1 non-relativistic Coulomb 
potential. Smaller systems (under N = 21 particles) indicate that RR state might be 
enhanced for cos 2 9 slightly less than 1 in (tt, 2) LL, which corresponds to a small relative 
decrease of V\ pseudopotential [8]. However, the largest system one can obtain by exact 
diagonalization, N = 21, suggests that the maximum overlap is practically at the bare 
Coulomb point. It cannot be determined with certainty whether this is an intrinsic 
feature of the RR state or a finite size effect that plagues N = 21 system. Finite-size 
effects are generally very strong at this filling factor and further theoretical studies, 
possibly employing density-matrix renormalization group techniques, are needed. On 
either side of the bare Coulomb point, the RR state is quickly destroyed. On the 
left, for cos 2 9 < 0.8 in (tt, 2) LL, there is an incompressible-compressible transition to 
a stripe phase. As mentioned earlier, states with broken translational symmetry are 
not adequately represented in the spherical geometry; however we can still detect the 
stripe by the presence of a sharp peak in the structure factor. One arrives at the same 
conclusion (but more transparently) in the torus geometry, where it is found that the 
system develops a manifold of nearly degenerate ground states for cos 2 9 < 0.8 that 
are characterized by a linear (uni-directional) array of pseudomomenta that defines the 
preferred direction of a stripe. As we move away from the pure Coulomb point in (7r, 1) 
LL, the ground state undergoes a change in shift for cos 2 9 > 0.15 in order to minimize 
the energy. The new value of the shift is that of a hierarchy state. Therefore, in chiral 
2DES, one can in principle study a subtle phase transition between two topological 
phases with Abelian and non-Abelian excitations. 
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Figure 4. (Color online) Charge gap (in units of e 2 /e^s) at v = 3/5 for the generalized 
material described by Eq. [U The gap smoothly decreases from the hierarchy phase 
down to zero (stripe phase). RR state appears in the transition region. The gaps are 
calculated for N = 18 particles at the shift —3. 



An interesting feature, shown in Fig. El is the presence of k = 3 RR correlations in 
a (2-7T, 1) LL. Note that the effective interaction in this case interpolates between pure 
n = and n = 2 (non-relativistic) LL Coulomb repulsion; therefore it is quite different 
from the n = 1 LL form-factor (FigfS]). In this case, RR state appears in a narrow 
window, between the stripe phase (for cos 2 ^ < 0.5) and the hierarchy state (shaded 
region in Fig. [31 for cos 2 # > 0.2). We note that N = 18 system has a significantly 
poorer overlap than all the other systems considered, including a larger N = 21 system, 
which leads us to conclude that the poor overlap might be a finite-size effecJJ More 
importantly, the excitation gap in the window where the ground state is well-described 
by the RR wavefunction is very small (and also smaller than in (rr, 1),(tt, 2) LLs). This, 
together with the presence of a strong incompressible hierarchy phase, would make it 
difficult to observe experimentally. 

In the remainder of this paper, we investigate the dependence of the excitation gap 
as we tune the parameters of the system. As mentioned above, the experimental gap 
at v = 12/5 is tiny It is therefore essential to find ways to enhance this gap in order 
to manipulate the state, e.g. via surface probes. To this end, we consider a generalized 
model for the material that involves a superposition of n = 0, n = 1 and n = 2 LL 
non-relativistic form factors [21] : 

F n {q) = cos 2 6L H (^^j+sm 2 9 cos 2 (f>L H+1 (^^j +sin 2 9 sin 2 0L w+2 ^ 

In Fig. H]we plot the charge gap at v — 3/5 for iV = 18 particles. The ground state 
of this particular system has a somewhat poorer overlap with the RR state than other 

% Such effects on the sphere often originate from the aliasing problem, but in this case it is not clear 
what the competing state is. Aliasing problem on the sphere was introduced in N. d'Ambrumenil and 
R. Morf, Phys. Rev. B 40, 6108 (1989). 
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Figure 5. (Color online) Charge gap (in units of e 2 /ei-^s) at v = 3/5 as a function of 
the screening parameter a and the distance of the dielectric, d/ls- The gap is largest 
for a close to unity and intermediate d ~ £g- For negative a, the gap vanishes and 
the system is in the stripe phase. Data is shown for N = 18 particles on the sphere at 
the shift —3. 

system sizes, however the charge gap shows no such difference from other cases. Along 
certain lines, indicated by arrows, the phase diagram reduces to one of the cases studied 
above. Note that proper finite-size scaling needs to be performed in order to get the 
correct values for the gap; however, it was previously found for filling factors v = 1/3 
and v = 1/2 [25j E3] that this rigorous analysis produces values that are roughly in 
agreement with the ones shown in Fig. HI At v — 3/5 the number of available system 
sizes is too small to perform such an extrapolation reliably, and the plot in Fig. H] should 
only be viewed qualitatively. We find a smooth transition between the hierarchy state 
and a stripe, with the gaps gradually dropping to zero for small cos 2 9. In the curved, 
narrow region around the transition, we previously found large overlaps with the RR 
state [23]. The excitation gaps in the transition region are much smaller than in the 
hierarchy phase, but they are nonetheless higher than at pure n — 1 non-relativist ic 
LL point (cos 2 9 = 0, cos 2 = 1). Therefore, tuning the band structure does allow for 
a modest enhancement of the gap in the RR phase. Qualitatively similar results are 
obtained for the neutral, instead of charge, gap. 

A larger variation of the gap can be achieved in the setup that utilizes the dielectric 
screening (Eq. |3]). In Fig. [5] we show the charge gap as a function of the screening 
parameter a and the distance of the dielectric plate, d/£ B . Unless a is close to unity, 
the gap does not vary significantly with d/is, which hints at the strong finite-size effects 
present even in this, fairly large, system of N = 18 particles. However, as a becomes 
close to unity, the gap displays a clear maximum for d ~ Eg, which is roughly two times 
larger than for the bare interaction. The overlap with the RR wavefunction is zero for 
a < —0.5, but fairly large otherwise. Therefore, at least for positive a, the gap can be 
tuned without reducing the good overlap with the model wavefunction. 
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4. Conclusion 

In this paper we have explored the stability of k = 3 RR state, the prime candidate 
for topological quantum computation, in the phase diagram as we tune the effective 
interaction via band structure parameters (that depend on the chirality of the material), 
or via dielectric screening. Numerical calculations suggest that the RR state is most 
stable in the vicinity of the bare n = 1 non-relativistic LL Coulomb potential, the 
same one which describes GaAs materials. Nevertheless, the chiral 2DESs offer several 
advantages over GaAs. Firstly, we find that RR state can be realized in various (as 
opposed to a single) LLs, where the interactions have different form from the one in 
GaAs. More importantly, the excitation gaps can be enhanced and the nature of the 
RR state can be further probed by driving in situ quantum phase transitions, either to 
compressible states with stripe ordering, or to an incompressible (hierarchy/composite 
fermion) state. 

5. Acknowledgements 

This work was supported by DOE grant DE-SC0002140. Y.B. was supported by the 
State of Florida. We acknowledge previous collaboration with R. Thomale. 

References 

[1] D. C. Tsui, H. L. Stormer, and A. C. Gossard, Phys. Rev. Lett. 48, 1559 (1982). 

[2] R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

[3] G. Moore and N. Read, Nucl. Phys. B 360, 362 (1991). 

[4] A. Yu. Kitaev, Ann. of Phys. 303, 2 (2003); ibid 321, 2 (2006); M. Freedman et at, Comm. Math. 

Phys. 227, 605 (2002); ibid 228, 177 (2002); C. Nayak et al, Rev. Mod. Phys. 80, 1083 (2008). 
[5] The Quantum Hall Effect, 2nd ed., edited by R. E. Prange and S. M. Girvin, Springer- Verlag, New 

York, 1990. 

[6] F. D. M. Haldane, Phys. Rev. Lett. 51, 605 (1983). 

[7] J. K. Jain, Composite fermions, (Cambridge University Press, 2007). 

[8] N. Read and E. Rezayi, Phys. Rev. B 59, 8084 (1999). 

[9] M. Levin, B. I. Halperin and B. Rosenow, Phys. Rev. Lett. 99, 236806 (2007); S.-S. Lee, S. Ryu, 
C. Nayak and M. P. A. Fisher, Phys. Rev. Lett. 99, 236807 (2007). 
[10] R. Willett et at, Phys. Rev. Lett. 59, 1776 (1987). 

[11] L. Hormozi, N. E. Bonesteel, and S. H. Simon, Phys. Rev. Lett. 103, 160501 (2009). 
[12] J. S. Xia et al, Phys. Rev. Lett. 93, 176809 (2004). 

[13] A. Kumar, G. A. Csathy, M. J. Mantra, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 105, 
246808 (2010). 

[14] P. Bonderson and J. K. Slingerland, Phys. Rev. B 78, 125323 (2008); Parsa Bonderson, Adrian E. 
Feiguin, G. Mooller, J. K. Slingerland. larXiv:0901.4965| A. Wojs, Phys. Rev. B 80, 041104(R) 
(2009). 

[15] E. H. Rezayi and N. Read, Phys. Rev. B 79, 075306 (2009). 

[16] Z. Papic, N. Regnault, and S. Das Sarma, Phys. Rev. B 80, 201303 (2009). 

[17] Michael R. Peterson, Kwon Park, and S. Das Sarma, Phys. Rev. Lett. 101, 156803 (2008). 

[18] R. H. Morf, Phys. Rev. Lett. 80, 1505 (1998). 

[19] E. H. Rezayi and F. D. M. Haldane, Phys. Rev. Lett. 84, 4685 (2000). 
[20] F. C. Zhang and S. Das Sarma, Phys. Rev. B 33, 2903 (1986). 



12 



[21] A. H. Castro Neto et ai, Rev. Mod. Phys. 81, 109 (2009). 

[22] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010). 

[23] B. A. Volkov and O. A. Pankratov, Pisma Zh. Eksp. Teor. Fiz. 42, 145 (1985); JETP Lett. 42, 

178 (1985); B. Buttner et ai, Nature Phys. 7, 418 (2011). 
[24] Z. Papic, D. A. Abanin, Y. Barlas, and R. N. Bhatt, Phys. Rev. B 84, 241306 (2011). 
[25] Z. Papic, R. Thomale, and D. A. Abanin, Phys. Rev. Lett. 107, 176602 (2011). 
[26] A. A. Koulakov, M. M. Fogler, and B. I. Shklovskii, Phys. Rev. Lett. 76, 499 (1996); R. Moessner 

and J. T. Chalker, Phys. Rev. B 54, 5006 (1996). 
[27] X. Du et ai, Nature 462, 192 (2009); K. Bolotin et ai, Nature 462, 196 (2009); C. R. Dean et ai, 

Nat. Phys. 7, 693 (2011). 
[28] K. S. Novoselov et ai, Nature 438, 197200 (2005); Y. Zhang et at, Nature 438, 201204 (2005). 
[29] K. Novoselov et ai, Nat. Phys. 2, 177 - 180 (2006). 

[30] V. M. Apalkov and T. Chakraborty, Phys. Rev. Lett. 97, 126801 (2006); C. Toke and J. Jain, 
Phys. Rev. B 75, 245440 (2007); Z. Papic, M. O. Goerbig, and N. Regnault, Phys. Rev. Lett. 
105, 176802 (2010). 

[31] W. Bao et ai, Nat. Phys. 7, 948 (2011); T. Taychatanapat et at, Nat. Phys. 7, 621 (2011); A. 

Kumar et ai, Phys. Rev. Lett. 107, 126806 (2011); C. H. Lui et ai, Nat. Phys. 7, 944 (2011). 
[32] Y. Barlas, K. Yang, and A. H. MacDonald, larXiv:1110.1069l 

[33] K. Nomura and A. H. MacDonald, Phys. Rev. Lett. 96, 256602 (2006); M. O. Goerbig, R. Moessner, 

and B. Doucot, Phys. Rev. B 74, 161407 (2006). 
[34] B. E. Feldman, J. Martin, A. Yacoby, Nature Phys. 5, 889 (2009). 

[35] Steven H. Simon, E. H. Rezayi, N. R. Cooper, and I. Berdnikov, Phys. Rev. B 75, 075317 (2007). 
[36] D. Yoshioka, B. I. Halperin, and P. A. Lee, Phys. Rev. Lett. 50, 1219 (1983); F. D. M. Haldane 
in Ref. [5]. 

[37] X. G. Wen and A. Zee, Phys. Rev. Lett. 69, 953 (1992). 

[38] M. Hermanns, Phys. Rev. Lett. 104, 056803 (2010); G. J. Sreejith, C. Toke, A. Wojs, and J. K. 

Jain, Phys. Rev. Lett. 107, 086806 (2011). 
[39] S. He, S. H. Simon and B. I. Halperin, Phys. Rev. B 50, 1823 (1994). 

[40] E. H. Rezayi, F. D. M. Haldane, and K. Yang, Phys. Rev. Lett. 83, 1219 (1999); F. D. M. Haldane, 

E. H. Rezayi, and K. Yang, Phys. Rev. Lett. 85, 5396 (2000). 
[41] B. Andrei Bernevig and F. D. M. Haldane, Phys. Rev. Lett. 100, 246802 (2008). 



